Dilute Hard "Sphere" Bose Gas in Dimensions 2, 4 and 5 
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The ground state energy for a dilute hard "sphere" Bose gas in various dimensions is studied 
theoretically. 

PACS numbers: 03.75.Hh, 05.30.Jp, 67.85. Be, 03.65.-w 



Consider a collection of N Bose hard spheres in a pe- 
riodic L X L X L box. The ground state energy Eq of the 
system was given [l| in an asymptotic expansion in 1957: 



-j^ = 47rap[l 



128 



(1) 



in the limit that iV — > oo at fixed density p — N/ft, for 
small value of pa^. We follow the notation of reference 
1, now usually called LHY. At that time equation (1) 
was a purely theoretical result. Today with the amazing 
developments of laser technology and computer power it 
has become possible [21 to test the validity of the expan- 
sion (1). In the present paper we attempt to find similar 
expansions for the same problem in dimensions 2, 4 and 
5. 

For dimension 1, the problem had been solved in ref- 
erence 3. 



REVIEW OF DIMENSION 3 

The one dimensional problem is special since there 
is no diffraction, only reflection, in one dimension. In 
higher dimensions diffraction is present. 

For the ground state only S-wave scattering need be 
y considered in the limit of fixed p, N ^ oo and a ^ 0. 
In this section, we review the derivation of equation (1) 
in dimension 3. 

The boundary condition that the wave function van- 
ishes when |ri — Tjl = is mathematically simple to de- 
fine, but difficult to analyze. We shall therefore replace it 
by a potential energy called the "pseudopotentiaF origi- 
nally [5| due to Fermi: 
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All equations with a square bracket refer to equations in 
LHY. 

To generalize to other dimensions we need to under- 
stand why the pseudopotential can replace the boundary 
condition when two spheres touch, at least in low orders 



of a. Now it is known from electrostatics that 



r 



Thus for 2 bodies 



if 



[-V^ + inaS^{r)—r]ip = ATraS'^ir)[—r - 1], 
or or 



ip = l 



(2) 



(3) 



(4) 



Now (a) the RHS of (3) when operating on any function 
that is not singular at r = gives zero, (b) (p satisfies the 
boundary condition that <y9 = at r = a, so that ip is the 
correct S wave scattered wave function for the collision 
of particles 1 and 2. This means that the N — 2 problem 
can be replaced by a pseudopotential 



c d 

V = 8'KaS-'{ri2)^ — m 
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(PP3) 



where we have taken into account the reduced mass for 
the 2 body system. For the N body problem we thus 
obtain the potential V of equation [2] in LHY. 

The next step is to use the simpler potential V' [3] in 
LHY which operating on any nonsingular wave function 
gives the same result as V. 

With V' LHY obtained for the ground state the first 
order perturbation energy and wave function. With V' 
the second order perturbation energy, however, diverges. 
LHY showed that that is due to the fact V is not V in 
the second order calculation. Using V and not V in the 
second order calculation, LHY showed that 



Eo = 47rap7V-^'[fc2+fco'-fc^F + 2fc§-^], [23] 

which is convergent and gave equation (1) above. 



DIMENSION 4 

To generalize to dimension 4 we need generalizations 
of equations (2)^(4): 



1 



-47r'^r(r). 



(5) 



V^ +A7r' a^S\v)^r^]^ = An^'S^r) 
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and 



= 1 - — 
The pseudopotential now becomes 



- 1], (6) 



(7) 



V = 8TT^a^S\r) 






(PP4) 



and 



V 



' = AtP-o^ /d'^rV'*(r)V'*(r)V'(r)V'(r) (8) 

We can now follow the steps in LHY to arrive at the 
energy E^: 

Eo^in'a'pN-J2'[k' + k',-k^k^ + 2^-^], (9) 

where now 

ko = v^STT^pa^. (10) 

The summation in (9) converts into an integral: 

NGAn^a^p' fd^[l+e-^Ve + 2-—], (11) 

Jo ^s 

where k — fcg^. This last integral has a ultra violet diver- 
gence: 



l^=-b<\°" = -bnkr- (12) 

2 ^ 2 ' 2 ' ^ ^ 

The cutoff for large k is ~ - . Thus the integral is equal 
to ' " 



Mn-r^. We have thus 

2 k-oa 

— = in^a^p + 32n*a^p^\lnkoa\ + 



(13) 



The depletion of the k = state for dimension 3 was 
given in LHY: 



(nk=o) =N[1- 



iV^ 



' pa-^ 



[40b] 



One arrives at a similar expression for dimension 4: 

(nk=o) =N[l-l3pa^\lnpa'\ + ---]. (14) 

where the logarithmic factor has the same origin as in 
(13) above. Here /3 is a numerical coefficient. Comparing 
equations (1) and [40b] for dimension 3 with equations 
(13) and (14) for dimension 4, we see that the expansion 
parameter changes from dimensions 3 to 4 as follows: 



' pa^ —f pa^\lnpa'^\. 



(15) 



DIMENSION 5 

Generalization of these perturbation calculations to di- 
mension 5 proceeds seemingly without problem, but what 
replaces (14) becomes 



(nk=o) = N[l - oo] 



(16) 



where cx) represents a linear divergence at large k. This 
means that the single particle state k = will be totally 
depleted by converting into many pairs (k, -k) of single 
particle states. In other words for interacting systems 
the BEC in dimensions 5 has [6| basic group of 2 instead 
of 1, unlike the BEC of dimensions 3 and 4. The details 
of this BEC need further study. 



DIMENSION 2 

In dimension 2 what replace (2) and (4) for dimension 
3 are 



and 



V^(-/nr) = -2tt5^{v). 



(p — In—, 
a 



(17) 



(18) 



if is the correct scattering amplitude, but it is impos- 
sible to generalize (3) to dimension 2 because unlike in 
dimensions 3 and 4, (18) does not approach a limit as 
r — > oo. However it changes very slowly at large r. To 
approximately solve this problem we search for a V which 
satisfies approximately — V^i/j = —Vip. Using (17) and 
(18) this equation become 



-2tt5^{v) 



-V[ln{r/a)]. 



(19) 



Now for a dilute system, the particles are at a large dis- 
tance of order 



R = Vp"'. 



(20) 



from each other. Thus we may put r = R in the RHS of 
(19) and get 

-2T:5^{v)'^-V[ln{R/a)]. 



or 
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(PP2) 



Here we take advantage of the fact that the logarithm 
varies very slowly with its argument if the latter is large. 
(The above derivation is not rigorous but I believe it to 
be correct.) 

With (PP2) we proceed as in LHY and obtain 



Eo 



AnpN 
\lnpa'^\ 



Y,'[k' + kl^k^k-^ + 2kl] (21 



Dimension 


Pseudopotential for fixed scatterer 


Eo/N 


Comments 


1 


— 


7rV/3(l-pa)^ 


NoBEC 


2 


^A-kS-' 


\ln{pa^)\ + 


Almost BEC 


\ln(pa^)\ 


3 


STraS^-§^r 


4.pa(l + i^ + ...) 


BEC. BG = f Boson 


4 


^-'a^^'l^r- 


47rVa^(l + 47rVa'*Kn (pa*)! + • ■ ■ ) 


BEC. BG = f Boson 


5 




? 


BEC. BG = pair of Bosons 



TABLE f: Comparison of BEC in dilute interacting Boson gases for different dimensions. (BG = basic group [6]) 



where 



fcg = %TTp/\lnpa 



There is no subtraction term in (21), unlike [23], be- 
cause now V' = V. The summation gives a contribution 

to Eq/N: 



Idirp 



InA 



where A — \lnpa'^\. But it is not clear that this term is 
reliable. 

We summarize the results for dimension n in Table I 
for dilute interacting Boson systems. 

With incredible new technology it has been possible to 
study experimentally 2D and ID systems [71] . Perhaps the 
qualitative and quantitative features exhibited in Table 
I above for different dimensions could be studied with 
these new technologies and with powerful computers. 
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